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Recent studies of the thermodynamic phase diagrams of the Gross-Neveu model (GN2), 
and its chiral cousin, the NJL2 model, have shown that there are phases with inhomogc- 
neous crystalline condensates. These (static) condensates can be found analytically because 
the relevant Hartree-Fock and gap equations can be reduced to the nonlinear Schrodinger 
equation, whose deformations are governed by the mKdV and AKNS integrable hierarchies, 
respectively. Recently, Thies et al have shown that time-dependent Hartree-Fock solutions 
describing baryon scattering in the massless GN2 model satisfy the Sinh-Gordon equation, 
and can be mapped directly to classical string solutions in A0IS3. Here we propose a geomet- 
ric perspective for this result, based on the generalized Weierstrass spinor representation for 
the embedding of 2d surfaces into 3d spaces, which explains why these well-known integrable 
systems underlie these various Gross-Neveu gap equations, and why there should be a con- 
nection to classical string theory solutions. This geometric viewpoint may be useful for higher 
dimensional models, where the relevant integrable hierarchies include the Davey-Stewartson 
and Novikov-Veselov systems. 

PACS numbers: ll.10.Kk, 11.15.Kc, ll.10.Lm, 02.40.Hw 

The Gross-Neveu model [IH3] is a remarkable 1 + 1 dimensional interacting fermionic model that 
shares some important features with quantum chromodynamics (QCD): it is asymptotically free, it 
has dynamical chiral symmetry breaking, and it has a limit of large number of flavors that behaves 
like the 't Hooft large N c limit of QCD. It has been well studied, but some surprising features 
have come to light only relatively recently. For example, at finite temperature and density, in the 
Nf — > 00 limit, the phase diagram shows regions in which the system prefers to form a spatially 
inhomogeneous crystalline condensate. This is true both of the Gross-Neveu model, GN2, that has 
a discrete chiral symmetry [H [5j , and of the N JL2 model that has a continuous chiral symmetry 
pi [7], although the phase diagrams are very different. The inhomogeneous phase of the GN2 
model has been verified on the lattice [8] . Similar Id inhomogeneous condensates have been found 
in related models, also in higher-dimensions [9H12|. The relevant gap equations can be solved 
analytically at finite temperature and density because of a rich integrability structure underlying 
the Gross-Neveu models. Important aspects of this integrability were recognized already in the zero 
temperature and zero density analyses [21 IT3HT5] , but the integrability properties have profound 



2 



implications for the analysis of these models at nonzero temperature and nonzero density \16\ [T7] . 
For example, the thermodynamic Ginzburg-Landau expansions of the GN2 and NJL2 models are 
expansions in functionals of the associated condensate, and these functionals are precisely the 
conserved quantities of the modified Korteweg-de Vries (mKdV) [IB] and Ablowitz-Kaup-Newell- 
Segur (AKNS) hierarchies [19], respectively; this permits the gap equation to be solved and re- 
summed to all orders. The (static) condensate satisfies the nonlinear Schrodinger equation, and the 
deformations of this equation encode the spectral information needed to study the thermodynamics; 
these deformations are governed by the mKdV and AKNS hierarchies, respectively, for real or 
complex condensates. 

More recently |20| , it has been shown that there are real time-dependent condensates that satisfy 
the Hartree-Fock equations of the GN2 system, and that these condensates are associated with the 
Sinh-Gordon equation, in light-cone coordinates. These solutions describe various baryon scattering 
processes [201 [21]; furthermore, they can be mapped directly [20] to various classical string solutions 
in AdS3 [22J. This is a much deeper result than that concerning static inhomogeneous condensate 
solutions, and so it raises the question of how this works in the gap equation language, and how 
this relates to the static condensates. Furthermore, at first sight it is quite surprising to see such 
an explicit connection between the GN2 model and string theory. While relations between string 
theory and Gross-Neveu models have been proposed previously [23ti25] , these have not been based 
on the inhomogeneous phases at finite temperature and density. 

In this paper we suggest a natural geometric explanation for the origin of the integrable Sinh- 
Gordon system in the Hartree-Fock analysis of the GN2 model, and our explanation makes it 
immediately clear why there is such a direct connection with classical AdS3 string solutions. Our 
proposal is strongly motivated by the seminal work of Regge and Lund [14J, who first explicitly 
recognized the geometric connection between string solutions, the sine-Gordon equation and the 
problem of embedding surfaces, and the work of Pohlmeyer [T3], and Neveu and Papanicolaou [15] . 
connecting sigma models, integrable systems and fermionic models. Our geometric explanation is 
based on the fact that the inhomogeneous condensates (static and time-dependent) of the GN2 
model all come from spinors that satisfy a nonlinear Dirac equation (in the terminology of [201 |2"T] 
they are "Type I": they solve the Hartree-Fock equations mode-by-mode, see Q below). This 
nonlinear Dirac equation has a simple geometrical interpretation in terms of the embedding of 
surfaces of constant mean curvature into 3 dimensional spaces. The key mathematical result is 
simple to state: it is a classical result in differential geometry [26-33J that the problem of immersing 
a 2d surface into a higher dimensional space may be expressed naturally in spinor language, with the 
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geometric Gauss-Codazzi equations appearing as the compatibility conditions for a Dirac equation. 
Initially suggested by Weierstrass and Enneper for minimal surfaces in IR 3 , this construction is now 
recognized as being far more general, even describing embeddings of surfaces into Lie groups and 
symmetric spaces. For surfaces of constant mean curvature, this Dirac equation becomes a nonlinear 
Dirac equation, and we show here that with appropriate choices of signature and manifolds, it is 
precisely the nonlinear Dirac equation of the GN2 model. Furthermore, the Sinh-Gordon equation 
appears because for constant mean curvature surfaces the square root of the induced metric factor 
of the emdedded surface is proportional to the spinor condensate, and the Gauss equation for this 
metric factor is precisely of Sinh-Gordon form. The relation to classical string solutions follows 
because constant mean curvature surfaces can be viewed as minimal surfaces [i.e. zero mean 
curvature], which are just string solutions, immersed in a space of constant curvature, such as 
AdS3. The relation to the static solutions found in (5j [TJ, [161 [TT] is that the dimensionally reduced 
version is the problem of immersing a curve [rather than a surface] into the 3d space, and this is 
also a classical problem in differential geometry, for which the familiar Frenet-Serret equations can 
be written in Dirac form, a fact that has been rediscovered many times in various physical and 
mathematical contexts, ranging from the motion of vortex filaments to optical solitons |34H40| . 

In Section I we review the Hartree-Fock and gap equation approaches to the Gross- Neveu model 
GN2. Section II describes the surface embedding problem for three different cases: (i) embedding 
a Euclidean signature 2d surface into M 3 ; (ii) embedding a Minkowski signature 2d surface into 
M 1,2 ; (iii) embedding a Minkowski signature 2d surface into AdS3. The first is the simplest case, 
which we review as an introduction, while we show that the Minkowski cases are the ones that are 
relevant for the Gross-Neveu model GN2. In Section III we describe the reduction to the previously 
known static inhomogeneous condensates. In Section IV we show how these ideas can be used to 
search for inhomogeneous condensates in the 2 + 1 dimensional Gross-Neveu model, a model for 
which much less is known. 

I. GROSS-NEVEU MODEL, HARTREE-FOCK AND GAP EQUATIONS 

The Gross-Neveu model (GN2) is a massless self-interacting fermionic theory in two dimensional 
spacetime with Lagrangian PQ: 

£ GN = i^>H + ~ W) 2 (1) 
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where g 2 is a coupling constant. This model is interesting because it is asymptotically free, it 
exhibits dynamical mass generation, and it has a (discrete) chiral symmetry that is spontaneously 
broken. We have suppressed the flavor indices as we work exclusively in the '"t Hooft limit" of 
Nj oo, with g 2 Nf fixed. The Dirac operator is <f) = j^d^, with Dirac matrices 7 and 7 , 
to be specified later (for the Euclidean and Minkowski cases). Note that we use the relativistic 
physics notation, ip = tp^ 7 . [A word of caution: the mathematical literature concerning the spinor 
representation of surface embeddings often uses 'tp to mean i/>' . The distinction is very important.] 
There are two equivalent methods to seek approximate solutions for the Gross-Neveu model. 
One method, the Hartree-Fock approach, is to seek self-consistent solutions to the Dirac equation 

(»0 -S)ip = (2) 

subject to the constraint 

(•0 tp) = ^2 ^v^v = — "2 S ( 3 ) 

states p ^ 

The simplest way [but not necessarily the only way, as we discuss later] to solve these equations is 
to find solutions ip p to the nonlinear Dirac equation for each mode p separately: 

[i$ -l$ p rp p ))ip p = (4) 

for some constant /, whose physical significance will become clear below. Thies et al refer to these 
mode-by-mode solutions to the nonlinear Dirac equation Q as "Type I" solutions |20|. I2T] . All 
known static and time-dependent solutions of the Hartree-Fock equations for GN2 are of this "Type 
I" form. The temperature (T) chemical potential (//) phase diagram of the massless GN2 model is 
characterized by a static crystalline condensate which has the form of a kink crystal: 

S(x) = m k sn(m x; k 2 ) (5) 

where sn is the usual Jacobi elliptic function, with elliptic parameter k 2 . With this S(x), the 
corresponding solutions ip p to the nonlinear Dirac equation can be constructed analytically, and 
the Hartree-Fock equations ([2]) and Q are satisfied [U [5]. This kink crystal condensate ^ is 
characterized by two parameters, m and k 2 , which together determine the amplitude, period and 
shape of the crystalline condensate. Thermodynamic minimization of the associated grand poten- 
tial determines m and k 2 as functions of T and /j,, so that the form of the crystal changes depending 
on temperature and density [H El E] • Indeed in one region of the phase diagram a vanishing con- 
densate is thermodynamically preferred (m = 0), while in another region a constant homogeneous 
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condensate is preferred (m ^ 0, and k 2 = 1), while in the crystalline phase the preferred condensate 
has and k 2 < 1. 

There are also exact time-dependent "Type I" solutions S(x, t) of the time-dependent Hartree- 
Fock equations, Q and Q, that also satisfy the nonlinear Dirac equation Q mode- by- mode [20] . 
The simplest nontrivial example is the boosted kink 

S(x, t) = tanh(2(a; - vt)/y/l - v 2 ) (6) 

from which one can construct a real kink-antikink scattering solution: 



v cosh(2x/\/l — v 2 ) — cosh(2f tj\J\ — v 2 ) 

t) = = — ( ' ) 

v cosh(2x/vl — v 2 ) + cosh(2ft/vl — v 2 ) 
The type I breather solution is not real |20j . Using the relation to the Sinh-Gordon equation, 
one can use Backlund transformations to construct more complicated solutions describing multi- 
baryon scattering processes |2LHI21|, Again, these time-dependent S(x, t) satisfy the time-dependent 
Hartree-Fock problem, and the corresponding solutions ip p to ^ can be constructed analytically, 
and satisfy the nonlinear Dirac equation Q for each mode p. 

A second, equivalent, method to approach the GN2 model is via a functional gap equation, 
where one replaces the fermion bilinear tpi/; by a bosonic condensate, S, leading to a gap equation 
(we suppress the coupling g): 

S(x,t) = j^—]ndet(i0-S(x,t)) (8) 



When S is constant, it is straightforward to solve this gap equation |41H43] . However, such a 
solution does not give the correct phase diagram [HOE] ■ When the condensate S is inhomogeneous, 
S = S(x,t), it is a highly nontrivial matter to solve the functional differential equation ([8]). For 
static, but spatially inhomogeneous, condensates S(x) it has been shown [16] that the gap equation 
can be mapped to the solution of the nonlinear Schodinger equation (NLSE), 

S" -2S 3 + cS = (9) 

where c is some constant. The most general bounded solution of ^ is given by the kink crystal 
(JHJ), with the identitifcation: c = m 2 {l + k 2 ). Thermodynamic minimization of the corresponding 
grand potential requires studying this solution as the two parameters m and k 2 are varied, which 
means we study deformations of the NLSE. These deformations are integrable and are governed 
by the mKdV hierarchy [HI Q2] , which means that the gap equation ^ is precisely the Novikov 
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equation for mKdV. More explicitly, in thermodynamic language, when S = S(x), we can expand 
the log det in a Ginzburg-Landau expansion: 



where a n (T,fx) are simple known functions of the temperature T and chemical potential [i, and 
a n [S] are known functionals of the condensate S(x) and its spatial derivatives [TBI [T7]. Indeed, 
the a n [S] are precisely the conserved quantities of the mKdV hierarchy. Using the integrability 
properties of the mKdV hierarchy, one finds that the solutions S(x) of the NLSE satisfy the gap 
equation Q term-by-term, in the sense that 



This is the gap equation analogue [16] of the mode-by-mode solution of the Hartree-Fock problem 
in @. 

In this gap equation language, the observation of Thies et al |20|, I21j . concerning the time- 
dependent solutions of the Hartree-Fock problem, is truly remarkable: it means that one has 
an analytic space- and time- dependent solution S(x,t) to the functional gap equation ([8]). At first 
sight, it is also surprising that these solutions should be directly related to classical string solutions. 
The connection to string solutions pointed out in [20\ [2Tj is via the Sinh-Gordon equation, which 
permits use of the explicit solutions constructed by Jevicki et al [22]. For these time-dependent 
condensates, S(x, t), the role of the NLSE for static solutions is played by the Sinh-Gordon equation. 
This has motivated us to look for the most natural language in which to understand these time- 
dependent solutions. We suggest here that the essence of the relation is due to an elementary 
property of differential geometry: the immersion of two dimensional surfaces is naturally expressed 
in terms of spinors and a Dirac-like equation, and for surfaces of constant mean curvature, this 
Dirac equation becomes a nonlinear Dirac equation. To match precisely to the GN2 model we have 
to choose carefully the surface and the space into which it is immersed, but the basic idea is due to 
this generalized Weierstrass representation of surface immersions, a classical result of differential 
geometry. 




(10) 



n 



Sa n [S] 
SS(x) 



oc S(x) 



for all n 



(11) 



II. NONLINEAR DIRAC EQUATION AND SURFACE EMBEDDINGS 



In this section we relate the nonlinear Dirac equation ([4]) to the problem of immersion of a 2d 
surface into a 3d space [261 [27]. We begin with a review of the Euclidean signature case, where the 
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most important result for our purposes is the generalized Weierstrass representation of a surface 
|28H33j . which states that the solution to the 2d (Euclidean) Dirac equation 

j(0 -5)^ = (12) 

defines a 2d surface in M 3 , such that the mean curvature H [defined below] of the surface is related 



to the potential S appearing in the Dirac equation (12) as: 



H = Wt (13) 



Thus, for a surface of constant mean curvature, H = I, the Dirac equation ( 12 ) becomes a nonlinear 
Dirac equation: 

i($ - I ^V) V 1 = (14) 

Notice that while this is a nonlinear Dirac equation, it is not of the Gross- Neveu form Q, where the 
nonlinear term is ipifi, instead of ip^ip. We then show that the nonlinear term ^jijj arises from em- 
bedding a Minkowski signature surface into a 3d space such as Mr> or AdS^. To establish notation 
and to compare with classical differential geometry, we first quickly review the standard Euclidean 
embedding case, the generalized Weierstrass representation, before turning to the Minkowski cases. 

A. Euclidean Signature: Conformal Immersion of 2d Surfaces in K 3 

1. Basic Differential Geometry of Euclidean Surface Embedding 

Consider a smooth, orientable surface X, parametrized by (x,y) € S C M 2 , embedded confor- 
mably into R 3 [261127] : 

X{x,y):E^R 3 (15) 

The conformal parametrization fixes the induced metric of the surface to be of the form ds 2 = 
f 2 (x, y)(dx 2 + dy 2 ), and in complex coordinates, z = x + iy, z = x — iy, we have 

ds 2 = f 2 (z,z)dzdz (16) 



The isometries of the metric (16) are holomorphic mappings z — > u(z). The surface can be 
described by the frame of three orthogonal vectors: the unit normal N to the surface, and the two 
derivatives, X x = d x X, and X y = d y X. It is more convenient to consider the complex tangent 



FIG. 1: Embedding of a 2d surface into K 3 . The surface can be characterized by the motion of the frame 
{N, X x ,X y }, consisting of the normal and tangent vectors, as it moves across the surface. This is described 
by the Gauss-Codazzi equations. 



vectors, X z = d z X, and X% = d z X , together with the the unit normal vector N. These three 



These inner products define the elements of the "first fundamental form". To characterize the 
surface, we consider how the frame changes as we move across the surface, which requires the 
second derivatives. In order to analyze the evolution of this moving frame, we expand X zz , X zz 
and X zz , in the basis {X Z ,X Z ,N}. The relevant basis projections are contained in the inner 
products 



which define the Hopf differential, Qdz 2 , and the mean curvature, H. The condition of commuta- 



vectors form a basis {e^} = {X Z ,X Z , N} in C 3 , with inner products: 

X Z .X Z — — , X Z .X Z — , X Z .X Z — 



(17) 



X ZZ .N = Q , X- Z - Z .N = Q* , X zz .N=-fH 



(18) 



tivity of the derivatives, euy 



zz 



%),zzi generates the Gauss-Codazzi equations: 



ffzz - fzfz + ~^H 2 f 4 




(19) 



Q 



z 



(20) 



Q* z 



(21) 



It is a simple exercise to show that these equations are invariant under isometries. 
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127] . The associated Gaussian curvature, K, which 



Functions f,Q and H satisfying the Gauss-Codazzi equations (19, 20 21) are necessary and 
sufficient to characterize a surface X in M 3 
is an intrinsic quantity, is given by 

4|Q| 2 



K = H Z 



-2f- 2 Hf 



(22) 



where in the last step we used the Gauss-Codazzi equation (19). Two important well-studied cases 



are surfaces of constant Gaussian curvature, and surfaces of constant mean curvature. In the first 



case, with _ff=constant, we see from (22) that the metric must satisfy the Liouville equation 



ln(/ 2 



-2Kf 



(23) 



In the second case, with H constant, H = I, we note that the Gauss-Codazzi equations (20, 21) 



imply that the Hopf differential is holomorphic, Q = Q(z), and therefore there exists an isometry 
where it is constant. Then, writing the metric factor as f 2 = 2 ^-e 9 , the first of the Gauss-Codazzi 



equations (19) becomes the Sinh-Gordon (ShG) equation: 



e 2S + 2l\Q\ sinh(#) = 



(24) 



2. Spinor Representation of Euclidean Surfaces 

The generalized Weierstrass representation expresses the Gauss-Codazzi equations in terms of 
a Dirac equation, and an associated zero-curvature condition. We choose the Euclidean gamma 
matrices, 70 = — iai, and 71 = — iui- Writing if) = (ipi,i'2) T , the Euclidean Dirac equation (12) 
becomes 



S -2d z 
2d- z s 



lf>l,z = ~\S 1p 2 
1p2,z = \Stpi 



(25) 



Note that if if) = (ifji,^) is a solution, and S is real, then (— V , 2?V'i) is a l so a solution. From 



these two solutions, form the 2x2 matrix <3?: 




(26) 



for which 



det <3? = if)^if) 



(27) 
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The connection to surface embeddings is expressed most compactly using the identification between 
vectors l£l 3 and quaternions 

X 3 X 1 - iX 2 



X = (Xt,X 2 ,X 3 ) o X 



Xi+iX 2 



-x. 



x a 



-tr (Xa a ) 



(28) 



(29) 



where a a are the Pauli matrices. Now, given the solution in (26), we construct the quaternionic 
representation of three vectors as 



X, 




X, 




N 



(30) 



where / = det<I>. It is straightforward to check that these define the moving frame of a surface, 
X = ^tr(Xff), and the geometric equations for the moving frame appear as the linear matrix 
equations 



$2 = V<5> 



(31) 



where U and V are 2x2 matrices. Using X zz = X zz we can express them in terms of the geometric 
quantities: 



U 



f 

H 



f 



V 



if o 



H 
2 

ii 
/ 



/ 



(32) 



By construction, the compatibility of the two linear equations in (31) implies that U and V satisfy 
the zero curvature condition, 



U- z -V z + [U, V] = 



(33) 



and this zero curvature condition encodes the Gauss-Codazzi equations (19, 20, 21). 



This construction is known as the generalized Weierstrass representation of the embedded sur- 



face, and all the geometrical information can be encoded in the Dirac equation (25 ). To summarize, 
we establish a dictionary between the Euclidean surfaces in IR 3 and spinors: 

Dirac equation : % (0 — S) ip = 
Mean curvature : H = S(^ip)~ l 
induced metric factor : / = det $ = t/j'tf) 



Hopf differential 



Q = O^UV'i -1P2 ^1,2) 

or = i?h^i-MU) 



(34) 



11 



Minimal surfaces (with H = 0) are described by a free massless Dirac equation, with S = 0, 
and the corresponding spinors are (anti)holomorphic. This is the original Weierstrass-Enneper 
representation. The more general constant mean curvature surfaces (with H = l) are described by 
a nonlinear Dirac equation, since H = I implies S = I ip^ip: 

Euclidean constant mean curvature surfaces «-)■ i{$ — I ip^ip)ip = (35) 



Furthermore, the Gauss equation (19) becomes the ShG equation 



Gauss equation -R- Sinh-Gordon equation : 9 ZZ + 2\Q\l sinh(#) = (36) 

where the metric factor f 2 is written as f 2 = e e . In terms of the function S(z, z) appearing in 
the Dirac equation, this ShG equation reads: 

SS z - z -S z S- z + S ^ = l 2 \Q\ 2 (37) 

Note that the general solutions of the ShG equation, and also the corresponding solutions of the 
Dirac equation, can be expressed in terms of theta functions 



B. Minkowski Signature: Immersion of 2d Surfaces in R 1 ' 2 
1. Basic Differential Geometry of Minkowski Surface Embedding 

To study the Gross-Neveu model in two-dimensional Minkowski spacetime it is more appropriate 
to use the light-cone coordinates, x± = x ± t, instead of the complex coordinates, z = x + iy, 
z = x — iy. So we now consider an embedding of the form X(x,t) : S C K 1 ' i— > R 1 ' 2 , with 
the conformal parameterization ds 2 = f 2 (—dt 2 + dx 2 ). The first fundamental form in light-cone 
coordinates is given by the inner products: 

X + .X_ = J — , X + .X + = , X_.X_=0 (38) 

To construct the Gauss-Codazzi equations, we define the Minkowski analogues of the Hopf differ- 
ential and the mean curvature, from the projections of the second derivatives: 

X ++ .N = -iQ {+) , X—.N = -iQ { - ) , X+-.N = - % -H f 2 (39) 
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The Gauss- Codazzi equations can easily be written by repeating the steps of the Euclidean con- 
struction in Minkowsi space: 



//+- - /+/- - \H 2 f 
(+) 



(+W-) 



Q 
Q 



(-) 



(40) 
(41) 
(42) 



As in the Euclidean case, surfaces of constant mean curvature {H = I) are characterized by a 
Sinh-Gordon equation. If H is constant, H = I, then we see that the Gauss-Codazzi equations 



(41, 42) imply that Q(±) are left/right moving. Furthermore, they can be set to be constants, 
by appropriate isometries. Defining f 2 



2y/\Q(- )Q (_>\ ^ e Q ^ er Q auss _c oc ; azz i equation (|40|) 



becomes the Sinh-Gordon (ShG) equation [note the opposite sign compared to (24)]: 



- 2l\ \QWQ(-)\ sinh(#) = 



(43) 



2. Spinor Representation of Minkowski Surfaces in 



(1,2 



The Minkowski version of the generalized Weierstrass representation follows analogously. But 
switching to Minkowski signature, the nonlinear Dirac equation of the Gross-Neveu model emerges 
in the constant mean curvature case. 

We define Minkowski Dirac matrices, 70 = 01, and 71 = ia2, and the Dirac equation is 



(i$ — S) ip = or 
which we can write as two equations: 



= 2 S ^ 2 



S -2id 4 



2id. 



V'2 



s 



tp = 



(44) 



(45) 



Now note that if ip = (^1,^2) is a solution of (44), and S is real, then ip = {—ip^^ip^F ls a l so a 



solution of (44). Thus, we now form the 2x2 matrix: 



>2 ^2 



(46) 



for which 



det $ = ipip 



(47) 
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These expressions should be contrasted with the Euclidean versions (26) and (27). The crucial 



difference from the Euclidean case is that the metric factor / will be identified with the condensate 
ijjij), instead of with the density tfi'ip. Ultimately this is due to the change of symmetry from SU(2) 
to SU(1, 1). The underlying symmetry is now 577(1, 1), and the basis vectors are parametrized by 
2x2 matrices together with the inner product XY = -\tx{XY): 



X + = if®' 1 | | $ 
1 0, 



X. 



++ 



N 




N = i$- L a 3 $ 



X+- ■ N 



(48) 



We use the same normalization condition, det <& = /. The motion of the basis frame is characterized 
by linear equations involving matrices U and V: 



(49) 



In terms of geometric quantities they have similar forms to the Euclidean case: 



U 




iff 



' f 




V 




(50) 



The Gauss-Codazzi equations (40, 41, 42) are encoded by the zero-curvature condition: 



U--V+ + [U, V] = 



(51) 



Before considering constant mean curvature surfaces, we give the geometric interpretation of 
the corresponding canonical energy momentum tensor of the Dirac system. In the t, x basis: 



1 (sh>- (Q {+) + Q { ~ ] ) ~(Q i+) - Q { ~ ] ) 

(q(+)_q(-)) W+(Q(+> + qH) 



(52) 



In terms of the lightcone coordinates, x±, the energy-momentum tensor components are related to 
the Hopf differentials, and Q^~\ and the mean curvature H: 

H 



Q 



(-) 



T 



(53) 



The energy momentum flow is identical with the Codazzi equations: 



(54) 
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Hence the Minkowski version of the surface/Gross- Neveu dictionary is: 

Dirac equation : (i$ — S)ip = 
mean curvature : H = S (^^j)^ 1 
induced metric factor : / = det <3? = iptp 

[q(+) = -^^ 1j+ -^, + ^i) , 
Hopf differentials : < (55) 

[Q(-)=t(^ v -^ 2 ) 

[In [20], the notation hi and /12 was used for these Hopf differentials: hi = —Q^ + \ and /12 = ■] 
Minimal surfaces (with H = 0) are described by a free massless Dirac equation, and the cor- 
responding spinors are right/left movers. Constant mean curvature surfaces (with H = I) are 
described by a nonlinear Dirac equation, since H = I implies S = I ipip: 

Minkowski constant mean curvature surfaces f-> (i$ — I '(Jnp)'ip = (56) 



Furthermore, the Gauss equation (40) becomes the ShG equation 



Gauss equation Sinh-Gordon equation : 6L| — 2\J Q( + ~)Q( ) I sinh(0) = (57) 



where the metric factor f 2 is written as f 2 = 2 VQ ( ^ e e _ j n ^ erms Q f ^j^g function S(x + ,x-) 
appearing in the Dirac equation, this ShG equation reads: 

S S + - - S+ S- - — = -Z 2 Q (+) Q (_) (58) 



Note the different signs relative to the Euclidean version (37). 



C. Minkowski Signature: Immersion of 2d Surfaces in AdS^ 

Our discussion so far gives a geometrical explanation of why the "Type I" solutions to the time- 
dependent Hartree-Fock problem must satisfy the Sinh-Gordon equation, as was argued in |20|. I21J . 
A further observation of \20\ [2Tj was that these solutions can be mapped directly to classical string 
solutions in AdS^. We now show that this follows from our geometrical argument because one 
can equivalently view the constant mean curvature surfaces in R 1 ' 2 as minimal surfaces (H = 0) 
in AdS^. Furthermore, classical string solutions in AdS^ are, by definition, minimal surfaces. 
For completeness, we give the details of this construction. We first describe the Gauss-Codazzi 
equations for a general conformal embedding M 1,1 1— > AdS^, and then the spinor formalism for 
minimal surfaces that correspond to classical string solutions. 
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1. Basic Differential Geometry of Minkowski Surface Embedding in AdS^ 

AdSs space can be can be realized as a hypersurface in M 2 ' 2 : 

X.X = -X 2 _ x -X$ + X\ + X\ = -R 2 (59) 
where R is the AdS radius. Therefore, the embedding M. > h-> AdS^ can be formulated naturally as 



an embedding M ' i— > R ' with the condition (59). In light-cone coordinates, where the induced 



metric is ds 2 = 2 Q (+) Q ( )R dx+dx^, the moving frame is constructed from the tangent vectors 
X± and two mutually orthogonal normal vectors, X and N. The elements of the first and second 
fundamental forms, and il//, are given by the inner products: 



x + .x_ = 2 Q (+) Q ( "^ 2 , x + .x + = x_.x_=o 



X ++ .N = -iQM , X ++ JV = -;qH , X + _.iV= 2 Q (+) y fl2g (60) 

Notice that we have chosen the metric factor differently than the previous cases. The reason for 
this choice is to make the mapping from minimal surface Gauss-Codazzi equations to the Lax 
equation explicit. With the same spirit of the previous sections, from the commutativity of second 
derivatives of the basis vectors, eu)^ = em |_ we obtain the Gauss-Codazzi equations: 

//+--/+/--£- f 2 = -Q i+) Q^R 2 (^ + h 2 ) 

qL +) = 2Q(+) Q ( ' )fl2 h + 

Q H = agwq^rf H _ (61) 

When the classical equations of motions 

AV - X± f^X = (62) 

of the string are imposed we have H = 0, since X^ is aligned in the direction of X. Therefore, a 

classical worldsheet of a string in AdS^ is a a minimal surface in the AdS-j space. Once again since 
Q(+) and QH 

are right/left moving, with an appropriate isometry we can set them to constants. 
With this choice, the first of the Gauss-Codazzi equations gives the ShG equation: 

//+--/+/_- ^ = -Q^Q^ = constant (63) 

which can be written as 
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where the metric factor is f 2 = 2\J Q( + )Q(~~) Re 9 . 

Comparing with the Mr' 1 i— > Mr' embeddings with constant mean curvature from the previous 
section, we observe a remarkable duality between the minimal surfaces embedded into AdS^, with 

Gaussian curvature K = 2/ 2 ln(/ 2 )_| , and constant mean curvature surfaces embedded into Mr' 2 , 

with Gaussian curvature K = —2f~ 2 ln(/ 2 )_| They are both characterized by the solutions of 

the Sinh-Gordon equation. Geometrically, it is a tradeoff between the constant intrinsic curvature 
i of AdSs, and the constant extrinsic curvature H. 



2. Spinor Representation of Minkowski Surfaces in AdS^ 

AdSs is associated with the group SO(2, 2), which can be factorized as SO(2, 2) = PSU(1, 1) x 
PSU(1, 1). It means that the embedding can be formulated by two spinors as opposed to one. 
This is consistent with the generalization of the Weierstrass representation for embeddings of 
surfaces into higher dimensional spaces, as discussed in |44t I45j . Following the notation of [20j, 
let us work with two complex coordinates, Z\ = X_i + zXq, and Z2 = X\ + 1X2, normalized as: 
— Z\Z\ + Z2Z2 = —R 2 . We introduce two sets of spinors, ip a and ipb, simultaneously normalized 
tpa'ipa = ^b^b = tyty ■ The spinor parameterization of the surface is given by: 

Z 1 = R^ , Z 2 = -R^M (65) 

Then the Gauss-Codazzi equations (or, equivalently, the classical string equations of motion plus 
the Virasoro constraints) translate into two Lax equations: 

lpa,+ = U a 1p a , tp a _ = V a lp a 

= U b yb b , ip b _ = V b ip b (66) 

They also fix the normalization to be ipip = f. 

- - S 

Ipa^a = ^b4>b = f = J (67) 

Here we identified the spectral parameter 1 with the AdS scale; I = ^. The Lax pair matrices with 
this fixed normalization are 

_ I ^ ^\ _( f 

S+ iQMl \ / q i/ 

Ub={ f J ■ Vb= (^ii 4 1 (68) 
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Consistency of the linear equations gives the zero-curvature conditions 

Ua,--Va, + + [U a ,V a ]=0 , U bt - - V bt+ + [U b , V b ] = (69) 



and these separately lead to the sinh-gordon equation (63). Physically, this Dirac system corre- 
sponds to the Gross- Neveu model since the "potential" is identified with the scalar interaction 
term S = I ipip. 

III. STATIC REDUCTION AND IMMERSIONS OF CURVES 

All the discussion so far has been searching for time-dependent solutions S(x,t) to the nonlinear 
Dirac equation Q. We have seen that S(x, t) must satisfy the ShG equation in the Gauss equation 
form [we consider here just the Minkowski case, with embedding into R ,2 ]: 

S S*+_ - S+ S- - ^S 4 = -I 2 Q {+) Q { - ] (70) 

To make connection with the previously found static solutions, S = S(x), to the gap equation, 
based on the NLSE Q, we note that in light-cone coordinates, assuming simple time dependence 
of the form e~ lEt , we can replace: 

d + ^^(d x -iE) , d.^±(d x + iE) (71) 



Thus, the Gauss equation (40) becomes 

SS"-(S') 2 -S 4 = -Al 2 QMQ^ (72) 

Now note that any solution of the NLSE ((oj), S" - 2S 3 + cS = 0, has a first integral (S") 2 = 
S 4 — cS 2 + d, where c and d are constants. For example, the kink crystal solution ^ satisfies the 
NLSE: S"-2S 3 + m 2 (l + k 2 )S = 0, and also (S') 2 = S 4 -m 2 (l + A; 2 )5 2 + m 4 A: 2 . Thus, any solution 
to the NLSE satisfies the ShG equation: 

S S" - OS") 2 - S 4 = -c (73) 



which is precisely the Gauss equation form of the ShG equation ( 72 ) . This explains the connection 
between the static condensates found previously in [T6j and the construction described here. It 
also explains the similarity between the generalized Weierstrass representation of the embedding 
of surfaces into M 3 and the Da Rios representation [3l] of the embedding of curves into M 3 , also 
known as the Hasimoto map [36J in the theory of the motion of vortex filaments. There is a vast 
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literature on the embedding of curves into 3d spaces, and it can be understood as a dimensional 
reduction of the surface embedding construction |30j. 

An important physical observation is that in studying the thermodynamics of the static con- 
densates, we actually need to know the corresponding fermionic energy spectrum, which translates 
geometrically into studying spectral deformations of the curves that preserve the form of the curve. 
It is known that the deformations of the NLSE are integrable and are governed by the mKdV in- 
tegrable hierarchy, and this explains why the thermodynamical Ginzburg-Landau expansion of the 
grand potential is an expansion in terms of the mKdV conserved quantities, as argued in |16} 117] . 



Indeed, if we re- write the Dirac equation (44) with the replacements (71), then we obtain a Dirac 
equation 

'* " ) ( | J | I (74) 

s -%d x J \ihj v*y 

which is precisely the Bogoliubov-De Gennes equation for the thermodynamics of the system with 



Dirac equation ( 44 ) . Geometrically, E is the deformation parameter generating the flow of isometric 



curves. In Section IV we apply this observation to the Gross- Neveu model in 2 + 1 dimensional 
spacetime. 

A. "Boosted" static potentials 

Given a static solution S = S(x), it is always possible to generate a class of time-dependent 
solutions S(x, t) by a Lorentz boost. Thus, in the case of Minkowski embedding, given a solution 



<S(x) to the static equation (73), with c rescaled to 1, we can define a new function 



S(x+, X-) = S ( a x + H — x\ (75) 

V a J 

which satisfies the (rescaled) time-dependent Sinh-Gordon equation 

SS+- - S + S- - S 4 = -1 (76) 

Writing a = e~ v , this is a boosted solution with velocity v = tanhr/: 

£(*+,*_) = S (^=(* ~vt)j (77) 

The spinor appropriately transforms as — > e azr, l 2 ^ . 

Similarly, for the Euclidean embedding into M 3 , the resulting one-dimensional equation, after 
rescaling, is 

SS" - {S'f + 5 4 = 1 (78) 
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Given a solution to this one- dimensional equation, we obtain a two-dimensional solution 

S(z,z) =S (az + ^z) (79) 
which satisfies the (rescaled) time-dependent Sinh-Gordon equation 

SS Z - Z -S Z S- Z + S A = 1 (80) 

Here the boost acts as a rotation. 

B. Examples 

At this section we recapitulate some known classical string solutions. It is useful to introduce 
the global coordinates (t, p, cf>): 

Zi = coshpe** Z 2 = sinhpe^ (81) 

We also label the worldsheet coordinates as (r, a) where the worldsheet metric is conformally 
-dr 2 + da 2 . 

1. Spinning string 

The classical solution of the Gubser-Klebanov-Polyakov string is given by 

p = p(a) t = KT (j) = lot (82) 

with the metric factor 



P 



/2 



4Q(+)g(-) 



= k 2 sn 2 (oja; k 2 ) = koj e~ 9 (83) 



S 2 

where ' means d a . Here the string is parameterized by the elliptic parameter k 2 that controls the 
length (hence the energy and angular momentum) of the string uj = 2K ^ \ and k = ku. The 
solution satisfies the one dimensional, static sinh-gordon equation: 

\6" = 6+- =uKs\nh9 (84) 

where AQ^Q^ = uj 2 k 2 . On the spinor side, it is identified with a fermion system subject to a 
periodic potential S = wns(wa; k 2 ). 
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2. Pulsating string 

The pulsating string 

p = p(r) t = kt (j) = ma (85) 

is described by the worldsheet with the metric factor S 2 = — (1 — /c 2 )w 2 nc 2 (u;T; k 2 ) and is associated 
with the one dimensional, time dependent sinh-gordon equation: 

-S + S 2 -S* = k 2 (l-k 2 )u 4 (86) 

IV. IMPLICATIONS FOR HIGHER-DIMENSIONAL GROSS-NEVEU MODELS 

In 2 + 1 dimensional spacetime the Gross-Neveu model is still renormalizable in the large Nf 
limit [46j . The phase diagram has been studied via a gap equation and via Hartree-Fock assuming a 
uniform condensate |47|, |4"8] . It has also been studied on the lattice assuming a uniform condensate 
|49j . It is still unclear whether there could be an inhomogeneous (static) crystalline condensate that 
has lower grand potential in some region of the temperature and chemical potential plane. This 
question was studied by Urlichs [38] , who argued that a stripe condensate of the form ^ satisfies 
the gap equation but is not thermodynamically preferred. So far, no truly 2d inhomogeneous 
condensate has been found that satisfies the gap equation. We argue here that if such a condensate 
is of "Type I" , then it should satisfy a ShG-type equation. 

The Lagrangian of the 2 + 1 dimensional Gross-Neveu (GN3) can be written 

a 2 

£ G N3 = + y (^ r 3 ^) 2 (87) 

(1 o\ 

where we choose the parity-preserving form of the interaction. Here, r 3 = I I , is a color 

v° -v 

matrix, so that \I/ is in fact a 4-component spinor. The associated Dirac equation is 

(»0 - A) * = (88) 

where the parity-even condensate is A = ^frs^. Since T3 is diagonal, we can choose to work with 
2-component spinors, ip, with both signs of a condensate S = tpip. A suitable choice of Dirac 
matrices in 2 + 1 dimensions is: 7 = 03, 7 1 = and j 2 = 102- To seek static solutions, we 
use the complex coordinates z = x + iy and z = x — iy for the spatial directions, so that the 
2-component Dirac equation becomes 

-2d, -S U 2 / W2/ 



21 



Again, the eigenvalue E will describe deformations, so to find the basic solution as an embedding of 
a surface into a 3d space, we first consider the E = equation. Now, observe that if if) = (ip\,4 ) 2) T 
is a solution, then so is ifi = (tp2,ipi) T , for S real. Thus, we form the 2x2 matrix <I> as before: 



S = (90) 




Note that 

det$ = (|^i| 2 - |V> 2 | 2 ) = i>il> (91) 

recalling that 7 = cr 3 in this representation. Thus, if we define the two dimensional static conden- 
sate S(z, z) = tpip, then we find that this condensate must satisfy the Sinh-Gordon equation: 

SS ZZ -S Z S Z -^S 4 = -\Q\ 2 (92) 



where Q = {tpi^z ~ ^2 i'lz)- Note that the signs in (92) are different from those appearing in 



the Euclidean embedding case ( 37 ) , but the same as those appearing in the Minkowski embedding 
case (58). 

This construction is the same form as for the embedding of a 2d surface into the 4d space M 2,2 
considered by Konopelchenko and Taimanov [444 EH] > where two copies of a 2 x 2 Dirac equation 
are needed. Physically, in our case of the GN3 model, these two copies correspond to the two copies 
making up the full 4-component spinor \P. So, for a constant mean curvature embedding of R 2 into 
M 2 ' 1 we have a nonlinear Dirac equation of precisely the correct form for a "Type I" solution to 
the GN3 gap equation. This static condensate S must satisfy the Gauss equation in its ShG form 



(92). Furthermore, to study the thermodynamics of such a solution we need spectral deformations 
of this equation, which are governed by the modified Novikov-Veselov (mNV) system [50j 151] . 
which generalizes the mKdV system that governed the thermodynamics of the static condensates 
in the 1 + 1 dimensional Gross-Neveu model |16| IT7] . The thermodynamic implications of this new 
geometrical perspective on the GN3 gap equation are currently being investigated. The step from 
Id to 2d is a large one, as is familiar from the intricacies of the inverse scattering method, but our 
result suggests that if progress can be made in the 2 + 1 dimensional system, then the modified 
Novikov-Veselov system, which is a special case of the Davey-Stewartson system, should play a 
significant role |50H52j . 
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V. CONCLUSIONS 



To conclude, we have shown that there is a natural geometric interpretation of the recent 
results of Thies et al [20\ [2"T] concerning time-dependent condensate solutions of the Hartree-Fock 
problem for the Gross-Neveu model GN2. The geometric interpretation is based on the fact that the 
problem of embedding a 2d surface of constant mean curvature into a 3d space can be expressed as 
a nonlinear Dirac equation [2UH33] ; and with the appropriate choice of signature of the surface and 
the embedding space, this nonlinear Dirac equation is precisely the nonlinear Dirac equation for 
mode-by-mode solutions of the GN2 Hartree-Fock equations. This perspective is strongly motivated 
by earlier work of Regge and Lund |14j , Pohlmeyer [13] , and Neveu and Papanicolaou |15j . For these 
Gross-Neveu solutions, the logarithm of the bilinear condensate automatically satisfies the Sinh- 
Gordon equation, which geometrically is the Gauss equation for the induced metric factor. The 
connection with classical string solutions observed in |20|, I2T] follows immediately once we identify 
a constant mean curvature solution in IR 1 ' 2 with a zero mean curvature surface (i.e., a classical 
string solution) embedded in a 3d space of constant negative curvature, AdS3. The relation to the 
explicit string solutions in [22] follows from the fundamental role of the Sinh-Gordon equation in the 
construction of string solutions |13fU5l[2"2"] l53 1 l54|. This geometrical perspective on the Gross-Neveu 
gap equation and Hartree-Fock equations permits us to address the currently unsolved question 
of the possibility of inhomogeneous condensates in the 2 + 1 dimensional Gross-Neveu model. We 
find that a static condensate that satisfies the Hartree-Fock equations in the same mode-by-mode 
manner, must satisfy a Sinh-Gordon equation. The thermodynamical implications of this result 
are currently being investigated. 

We have not yet addressed the so-called "Type II" solutions of the massive GN2 system or of the 
massless NJL2 system, each of which satisfies a generalized nonlinear Dirac equation. Assuming a 
nonlinear Dirac equation, with a mode-by-mode solution (i.e. a "Type I" solution), the NJL2 system 
leads to a clear geometrical interpretation |13H15| . But it is known that the static condensates for 
this system are not necessarily mode-by-mode solutions of the nonlinear Dirac equation, except 
for the simple chiral spiral solution A(x) = Ae l ^ x (which, incidentally, is the thermodynamically 
preferred solution of the massless NJL2 system). These more general condensates satisfy more 
complicated nonlinear equations, and it would be interesting to investigate their possible geometric 
form, along with the connection to brane- world scenarios of Gross-Neveu models [55J. We also note 
that nonlinear Dirac equations have been used recently in the analysis of Bose-Einstein condensates 
in a two-dimensional honeycomb lattice |56j. It would be interesting to see if this geometrical 
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approach could be useful in such systems. Another curious relation worth mentioning is the fact 
that the Gauss-Codazzi equations are also naturally related to the Hitchin equations |57U58j . which 
in turn arise in the relation between self-dual Chern-Simons solitons and integrable models [59Vl62j . 
This is a reflection of the sigma model form of the Pohlmeyer reduction problem [54] . 

VI. ACKNOWLEDGMENTS 

We thank the DOE for support under grant DE-FG02-92ER40716, and M. Thies for helpful 
correspondence. 



[1] D. J. Gross and A. Neveu, "Dynamical Symmetry Breaking In Asymptotically Free Field Theories," 
Phys. Rev. D 10, 3235 (1974). 

[2] R. F. Dashen, B. Hasslacher and A. Neveu, "Semiclassical Bound States In An Asymptotically Free 
Theory," Phys. Rev. D 12, 2443 (1975). 

[3] J. Feinberg, "All about the static fermion bags in the Gross-Neveu model," Annals Phys. 309, 166 
(2004) [arXiv:hep-th/0305240]. 

[4] M. Thies and K. Urlichs, "Revised phase diagram of the Gross-Neveu model," Phys. Rev. D 67, 
125015 (2003) [arXiv:hep-th/0302092]; M. Thies, "Analytical solution of the Gross-Neveu model at 
finite density," Phys. Rev. D 69, 067703 (2004) [arXiv:hep-th/0308164]. 

[5] M. Thies, "From relativistic quantum fields to condensed matter and back again: Updating the Gross- 
Neveu phase diagram," J. Phys. A 39, 12707 (2006) [arXiv:hep-th/0601049]. 

[6] V. Schon and M. Thies, "Emergence of Skyrme crystal in Gross-Neveu and 't Hooft models at finite 
density," Phys. Rev. D 62, 096002 (2000) [arXiv:hep-th/0003195]; V. Schon and M. Thies, "2D model 
field theories at finite temperature and density," arXiv:hep-th/0008175, in At the frontier of particle 
physics, vol. 3, page 1945, Festschrift in honor of Boris Ioffe, edited by M. Shifman (World Scientific, 
2000). 

[7] G. Basar, G. V. Dunne and M. Thies, "Inhomogeneous Condensates in the Thermodynamics of the 
Chiral NJL2 model," Phys. Rev. D 79, 105012 (2009) [arXiv:0903.1868 [hep-th]]. 

[8] P. de Forcrand and U. Wenger, "New baryon matter in the lattice Gross-Neveu model," PoS LAT2006, 
152 (2006) [arXiv:hep-lat/0610117]. 

[9] D. Nickel and M. Buballa, "Solitonic ground states in (color-) superconductivity," Phys. Rev. D 79, 
054009 (2009) [arXiv:0811.2400 [hep-ph]]; D. Nickel, "How many phases meet at the chiral critical 
point?," Phys. Rev. Lett. 103, 072301 (2009) [arXiv:0902.1778 [hep-ph]]; "Inhomogeneous phases in 
the Nambu-Jona-Lasino and quark-meson model," Phys. Rev. D 80, 074025 (2009) [arXiv:0906.5295 
[hep-ph]]. 



24 



[10] L. Y. Glozman and R. F. Wagenbrunn, "Chirally symmetric but confining dense and cold matter," 
Phys. Rev. D 77, 054027 (2008) [arXiv:0709.3080 [hep-ph]]; "Second order chiral restoration phase 
transition at low temperatures in quarkyonic matter," arXiv:0805.4799 [hep-ph]. 

[11] T. Kojo, Y. Hidaka, L. McLerran and R. D. Pisarski, "Quarkyonic chiral spirals," AIP Conf. Proc. 
1257, 732 (2010); T. Kojo, R. D. Pisarski and A. M. Tsvelik, "Covering the Fermi Surface with Patches 
of Quarkyonic Chiral Spirals," Phys. Rev. D 82, 074015 (2010) [arXiv: 1007.0248 [hep-ph]]. 

[12] I. E. Frolov, V. C. Zhukovsky and K. G. Klimcnko, "Chiral density waves in quark matter within 
the Nambu-Jona-Lasinio model in an external magnetic field," Phys. Rev. D 82, 076002 (2010) 
[arXiv: 1007.2984 [hep-ph]]. 

[13] K. Pohlmeyer, "Integrable Hamiltonian Systems And Interactions Through Quadratic Constraints," 
Commun. Math. Phys. 46, 207 (1976). 

[14] F. Lund and T. Regge, "Unified Approach To Strings And Vortices With Soliton Solutions," Phys. 
Rev. D 14, 1524 (1976); F. Lund, "Note On The Geometry Of The Nonlinear Sigma Model In Two- 
Dimensions," Phys. Rev. D 15, 1540 (1977); F. Lund, "Solitons and Geometry", In Nonlinear Equations 
In Physics and Mathematics, (Proceedings, Istanbul 1977), A. O. Barut (Ed), (D. Reidel, Boston, 1978). 

[15] A. Neveu and N. Papanicolaou, "Integrability Of The Classical Scalar And Symmetric Scalar- 
Pseudoscalar Contact Fermi Interactions In Two-Dimensions," Commun. Math. Phys. 58, 31 (1978). 

[16] G. Basar and G. V. Dunne, "Self-consistent crystalline condensate in chiral Gross-Neveu and 
Bogoliubov-de Gennes systems," Phys. Rev. Lett. 100, 200404 (2008) [arXiv:0803.1501 [hep-th]]; 
G. Basar and G. V. Dunne, "A Twisted Kink Crystal in the Chiral Gross-Neveu model," Phys. Rev. 
D 78, 065022 (2008) [arXiv:0806.2659 [hep-th]]. 

[17] F. Correa, G. V. Dunne and M. S. Plyushchay, "The Bogoliubov/de Gennes system, the AKNS 
hierarchy, and nonlinear quantum mechanical supersymmetry," Annals Phys. 324, 2522 (2009) 
[arXiv:0904.2768 [hep-th]]. 

[18] F. Gesztesy and H. Holden, Soliton Equations and their Algebro- Geometric Solutions, (Cambridge 
University Press, 2003) 

[19] M. J. Ablowitz, D. J. Kaup, A. C. Newell and H. Segur, "The Inverse scattering transform fourier 
analysis for nonlinear problems," Stud. Appl. Math. 53, 249 (1974). 

[20] A. Klotzek and M. Thies, "Kink dynamics, sinh-Gordon solitons and strings in AdS(3) from the Gross- 
Neveu model," J. Phys. A 43, 375401 (2010) [arXiv:1006.0324 [hep-th]]. 

[21] C. Fitzner and M. Thies, "Exact solution of an N baryon problem in the Gross-Neveu model," 
arXiv:1010.5322 [hep-th]. 

[22] A. Jevicki, K. Jin, C. Kalousios and A. Volovich, "Generating AdS String Solutions," JHEP 0803, 032 
(2008) [arXiv:0712.1193 [hep-th]]; A. Jevicki and K. Jin, "Solitons and AdS String Solutions," Int. J. 
Mod. Phys. A 23, 2289 (2008) [arXiv:0804.0412 [hep-th]]. 

[23] E. Antonyan, J. A. Harvey and D. Kutasov, "The Gross-Neveu model from string theory," Nucl. Phys. 
B 776, 93 (2007) [arXiv:hep-th/0608149]. 



25 



[24] A. Basu and A. Maharana, "Generalized Gross-Neveu models and chiral symmetry breaking from string 

theory," Phys. Rev. D 75, 065005 (2007) [arXiv:hep-th/0610087]. 
[25] J. L. Davis, M. Gutperle, P. Kraus and I. Sachs, "Stringy NJL and Gross-Neveu models at finite density 

and temperature," JHEP 0710, 049 (2007) [arXiv:0708.0589 [hcp-th]]. 
[26] L. P. Eisenhart, A Treatise on the Differential Geometry of Curves and Surfaces, (Ginn Co., Boston, 

1909). 

[27] H. Hopf, Differential Geometry in the Large, Lcct. Notes Math. 1000, (Springer, Berlin, 1983). 
[28] K. Kenmotsu, "Weierstrass formula for surfaces of prescribed mean curvature ", Math. Ann. 245, 89 
(1979). 

[29] A. I. Bobenko, "Integrable surfaces", Func. Anal. Appl. 24, 227 (1990); "Constant mean curvature 
surfaces and integrable equations", Russ. Math. Surv. 46, 1 (1991); "All constant mean curvature 
tori in R3, S3 and H3 in terms of thcta-functions" , Math. Ann. 290, 209 (1991); "Exploring Sur- 
faces through Methods from the Theory of Integrable Systems: Lectures on the Bonnet problem," 
[arXiv:math/9909003], in Surveys on Geometry and Integrable Systems, Advanced Studies in Pure 
Mathematics 51, 1 (2008). 

[30] B. G. Konopelchenko, "Induced surfaces and their integrable dynamics", Stud. Appl. Math. 96, 9 
(1996); B. G. Konopelchenko and I. A. Taimanov, "Constant mean curvature surfaces via integrable 
dynamical system", J. Phys. A 29, 1261 (1996) [arXiv:dg-ga/9505006]. 

[31] A.S. Fokas and I.M. Gel'fand, "Surfaces on Lie Groups, on Lie Algebras, and Their Integrability" , 
Commun. Math. Phys. 177, 203 (1996); A. S. Fokas, I. M. Gelfand, F. Finkel and Q. M. Liu, "A 
formula for constructing innitely many surfaces on Lie algebras and integrable equations", Sel. Math., 
New ser. 6, 347 (2000). 

[32] I. A. Taimanov, "Two-dimensional Dirac operator and the theory of surfaces", Russ. Math. Surv. 61, 
79 (2006). 

[33] S. P. Novikov and I. A. Taimanov, Modern geometric structures and fields, (AMS, Providence, 2006). 
[34] L. S. Da Rios, "Moto dun liquido indefinite con un filetto vorticoso di forma qualunque (On the motion 

of an unbounded liquid with a vortex filament of any shape)" , Rend. Circ. Mat. Palermo 22, 117 (1906); 

See R. L. Ricca, "Rediscovery of the Da Rios Equations", Nature 352, 561 (1991). 
[35] G. L. Lamb, "Analytical Descriptions of Ultrashort Optical Pulse Propagation in a Resonant Medium," 

Rev. Mod. Phys. 43, 99 (1971); G. L. Lamb, "Solitons On Moving Space Curves," J. Math. Phys. 18, 

1654 (1977); G. L. Lamb, "Solitons And The Motion Of Helical Curves," Phys. Rev. Lett. 37, 235 

(1976). 

[36] R. Hasimoto, "A soliton on a vortex filament", J. Fluid Mcch. 51, 477 (1972). 

[37] A. Sym and J. Corones, "Lie Group Explanation Of Geometric Interpretations Of Solitons," Phys. Rev. 

Lett. 42, 1099 (1979); A. Sym, "Soliton Surfaces", Lett. Nuovo Cim. 332, 394 (1982). 
[38] R. K. Dodd, "Soliton Immersions", Commun. Math. Phys. 197, 641 (1998). 

[39] A. Calini and T. Ivey, "Connecting geometry, topology and spectra for finite-gap NLS potentials", 



26 



Physica D 152, 9 (2001). 

[40] P. G. Grinevich and M. U. Schmidt, "Closed curves in R3 : a characterization in terms of curvature and 
torsion, the Hasimoto map and periodic solutions of the filament equation", Preprint: dg-ga/9703020. 

[41] U. Wolff, "The phase diagram of the innnite-N Gross-Neveu Model at finite temperature and chemical 
potential", Phys. Lett. 157B, 303 (1985). 

[42] T. F. Treml, "Dynamical Mass Generation In The Gross-Neveu Model At Finite Temperature And 
Density," Phys. Rev. D 39, 679 (1989). 

[43] A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini and R. Gatto, "Thermodynamics Of The Massive 
Gross-Neveu Model," Phys. Rev. D 51, 3042 (1995) [arXiv:hcp-th/9406117]. 

[44] B.G. Konopelchcnko and G. Landolfi, "Generalized Weierstrass representation for surfaces in mul- 
tidimensional Riemann spaces", Journ. Gcom. Phys. 29, 319 (1999) [arXiv:math/9804144]; B.G. 
Konopelchcnko, "Weierstrass representations for surfaces in 4D spaces and their integrable deformations 
via DS hierarchy", Ann. Global Anal. Geom. 16, 61 (2000) [arXiv:math/9807129]. 

[45] I. A. Taimanov, "Surfaces in the four-space and the Davey-Stewartson equations", J. Geometry and 
Physics 56, 1235 (2006), [arXiv:math/0401412]. 

[46] B. Rosenstein, B. Warr and S. H. Park, "Dynamical symmetry breaking in four Fermi interaction 
models," Phys. Rept. 205, 59 (1991); G. Gat, A. Kovner and B. Rosenstein, "Chiral phase transitions 
in d = 3 and renormalizability of four Fermi interactions," Nucl. Phys. B 385, 76 (1992). 

[47] K. G. Klimenko, "Three-dimensional Gross-Neveu model at nonzero temperature and in an external 
magnetic field," Z. Phys. C 54, 323 (1992). 

[48] K. Urlichs, 'Baryons and baryonic matter in four-fermon interaction models", PhD Thesis, Univ. Er- 
langen, 2007, and unpublished notes (2007). 

[49] S. Hands, S. Kim and J. B. Kogut, "The U(l) Gross-Neveu model at nonzero chemical potential," 
Nucl. Phys. B 442, 364 (1995) [arXiv:hep-lat/9501037] ; S. Hands, "Four fcrmion models at non-zero 
density," Nucl. Phys. A 642, 228 (1998) [arXiv:hep-lat/9806022]. 

[50] S. P. Novikov and A. P. Veselov, "Finite-zone two-dimensional periodic Shrodinger operators: potential 
operators", Dokl. Akad. Nauk SSSR 279, 784 (1984); "Two-dimensional Schrodinger operator: inverse 
scattering and evolutional equations", Physica D 18, 267 (1986); "Exactly solvable two-dimensional 
Schrodinger operators and Laplace transformations", Solitons, geometry, and topology: on the cross- 
road, Amer. Math. Soc. Transl. Ser. 2, 179, 109 (1997) [arXiv:math-ph/0003008vl]. 

[51] L. V. Bogdanov, "Veselov-Novikov equation as a natural two-dimensional generalization of the 
Kortcwcg-Dc Vrics equation", Theor. Math. Phys. 70, 219 (1987); L. V. Bogdanov, "On the two- 
dimensional Zakharov-Shabat Problem", Theor. Math. Phys. 72, 790 (1988). 

[52] I. A. Taimanov, "Modified Novikov- Veselov equation and differential geometry of surfaces", Amer. 
Math. Soc. Transl., Ser. 2, 179, 133 (1997) [arXiv:dg-ga/9511005]. 

[53] H. J. De Vega and N. G. Sanchez, "Exact Integrability Of Strings In D-Dimensional De Sitter Space- 
Time," Phys. Rev. D 47, 3394 (1993). 



27 



[54] M. Grigoriev and A. A. Tscytlin, "Pohlmcycr reduction of AdS$ x S 5 superstring sigma model," Nucl. 

Phys. B 800, 450 (2008) [arXiv:0711.0155 [hep-th]]. 
[55] W. Bietenholz, A. Gfeller and U. J. Wiese, "Dimensional reduction of fcrmions in brane worlds of the 

Gross-Neveu model," JHEP 0310, 018 (2003) [arXiv:hep-th/0309162]. 
[56] L. H. Haddad and L. D. Carr, "The Nonlinear Dirac Equation in Bose-Einstein Condensates: 

Relativistic Linear Stability Equations, Nonlinear Localized Modes, and Cherenkov Radiation", 

arXiv:1006.3893v2; "The Nonlinear Dirac Equation in Bose-Einstein Condensates: Foundation and 

Symmetries", Physica D 238, 1413 (2009) [arXiv:0803.3039vl]. 
[57] R. S. Ward, "Intcgrable and solvable systems, and relations among them," Phil. Trans. Roy. Soc. Lond. 

A 315, 451 (1985). 

[58] N. J. Hitchin, "The Selfduality Equations On A Riemann Surface," Proc. Lond. Math. Soc. 55, 59 
(1987). 

[59] B. Grossman, "Hierarchy of soliton solutions to the gauged nonlinear Schrodingcr equation on the 
plane," Phys. Rev. Lett. 65, 3230 (1990). 

[60] G. V. Dunne, R. Jackiw, S. Y. Pi and C. A. Trugenberger, "Selfdual Chern-Simons solitons and two- 
dimensional nonlinear equations," Phys. Rev. D 43, 1332 (1991) [Erratum-ibid. D 45, 3012 (1992)] 
[Phys. Rev. D 45, 3012 (1992)]. 

[61] G. V. Dunne, Selfdual Chern-Simons Theories, Lect. Notes Phys. M36, 1 (Springer, Heidelberg, 1995). 

[62] L. Martina, Kur. Myrzakul, R. Myrzakulov and G. Soliani, "Deformation of surfaces, integrable systems, 
and Chern-Simons theory", Journ. Math. Phys. 42, 1397 (2001). 



